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We present a theoretical analysis of the paradigm of encoded universality, using a Lie algebraic 
analysis to derive specific conditions under which physical interactions can provide universality. We 
discuss the significance of the tensor product structure in the quantum circuit model and use this to 
define the conjoining of encoded qudits. The construction of encoded gates between conjoined qudits 
is discussed in detail. We illustrate the general procedures with several examples from exchange- 
only quantum computation. In particular, we extend our earlier results showing universality with 
the isotropic exchange interaction to the derivation of encoded universality with the anisotropic 
exchange interaction, i.e., to the XY model. In this case the minimal encoding for universality is 
into qutrits rather than into qubits as was the case for isotropic (Heisenberg) exchange. We also 
address issues of fault-tolerance, leakage and correction of encoded qudits. 

I. INTRODUCTION 

Many machines have been invented to perform certain specific computational tasks. However in the language 
of computer science, these machines are not necessarily computers, meaning that they can not perform all possible 
computational tasks, i.e., they are not universal. What is considered as all possible computations depends on the 
underlying laws that govern the machines we use. In the classical world we can cast any computational task into a 
boolean function of the input (in bits). The evaluation of a computable function can then be reduced to a sequence 
of elementary logical operations acting on a few bits each: Classical circuits can be build from a universal gate set. 
The computational task of a quantum computer is to accurately perform any unitary operation on the input (qubits). 
These unitàries will have to be constructed from the bàsic units the computer will be equipped with, i.e., from a 
universal quantum gate set. Moreover these realistic quantum gates have to be physically realizable - they should 
involve only a few qubits at a time, for example. 

Research efforts in recent years has shown that realistic universal quantum gate sets exist and that uniform quantum 
circuits can be built. A physical realization needs to implement such a universal set in order to be considered a universal 
quantum computer. Some of the gates in such a set are easier to implement than others. In several proposals some 
two-qubit gates are much easier to implement than one-qubit gates, yet the latter ones are needed to complete a 
universal set of gates. We have previously presented results on how to make do with a restricted set of gates that 
is not universal. This restricted set does not accurately enact every unitary operation on the system qubits, yet 
can nevertheless provide universality on a restricted subspace [0J^]. This approach, that we have termed "encoded 
universality" , employs the philosophy of encoding the information, in order to make an a priori non-universal set of 
gates universal on the code ||| . This potentially reduces the experimental efforts of implementing gates that are hard 
to build. 

In this paper we show that our recent demonstration of universality from a single physical interaction achieved 
with the Heisenberg exchange interaction Q can be extended to an anisotropic interaction, namely to the XY model. 
This wcll-known Hamiltonian model is relevant to a number of two-dimensional condensed matter phenomena. The 
result shown here indicates the generality of the approach of seeking universal encodings tailored to specific physical 
interactions ||. 

A second aim of this paper is to expand the Lie algebraic analysis presented bricfly in Ref . j| , to show specifically the 
conditions under which selected elements of the àlgebra generated by the physical interactions can provide universality. 
We distinguish between the general results of universality on a d-dimensional subspace provided by this approach [|| , 
and a more practical attempt to develop encoded one- and two-qubit gates that implicitly imposes the Standard model 
of arbitrary one-qubit rotations and CNOT gates on the encoded qubits jï]|2| . The construction of these encoded gates 
is discussed in detail, including efficient implementation in terms of discrete operations jl| . We will also address the 
important issue of fault-tolerance and leakage detection and correction. This will furnish a complete picture of how 
fault-tolerance meshes with encoded universality. 

The structure of this paper is as follows. We shall first provide a brief summary of the concept of universality 
within the quantum circuit model of computation in Section [01. We then present the essentials of the Lie àlgebra 
analysis underlying the problem of simulation of a quantum circuit, and discuss the various levels of tensor product 
structure that are required to implement a quantum circuit structure. The representation theoretic framework is 
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presented in Section H], illustrated by a summary of the earlier result for isotropic exchange [§|J,mL We also provide 
the mathematical connection to results for universal computation on decoherence-free subsystems J2| . The core of the 
paper follows in Section IV, with a detailed presentation of the derivation of universal encodings for the anisotropic 
exchange interaction. Finally we address the issues of fault-tolerance and leakage (Sec. 0) and close with discussion 
of some open tòpics (Sec. VI). 



II. UNIVERS ALITY 



A. Quantum Circuits and Complexity 

To build a computer from a physical system for algorithmic applications means to provide its user with a set of 
"units" from which he or she can implement an algorithm. A classical computer has to evaluate boolean functions. A 
convenient model for classical computation is the uniform circuit-model. The input to the circuit is given by a string 
of bits; wires carry them to elementary gates (such as AND and NOT) which are again interconnected by wires. The 
computation carries the bits along the wires and executes the gates, when they are reached. The output wires carry 
the results. These circuits can be parametrized by the size of the problem, i.e. the number n of input bits. Informally 
speaking such a family of circuits is uniform if given n, the corrcsponding circuit can be constructed efhciently (i.e. 
in polynomial timc in n). This assures, that the complexity of a problem is not hidden into the design of the circuit. 
Such uniform circuits can be simulated at polynomial cost by a universal circuit, one of which is the circuit with the 
NAN D and the COPY gates only ||. The notion of complexity defines which problems are hard and which are 
easy according to the amount of space and time (i.e. gates) a circuit uses to solve them. Although classical universal 
circuits may differ in the bàsic gate-sets they use, any universal circuit can simulate any other efficiently, i.e. with 
only polynomial overhead in both space and time. Consequently the definition of complexity of a circuit can be made 
independently of the specifics of the gate-set used. 

The Quantum Circuit Model (QCM) can be viewed as a transposition of the classical circuit model to the quantum 
world. A quantum computer implcmcnts unitary operations. The input is given by a string of qubits, conncctcd 
to elementary gates. These gates constitute a universal set if the circuit can implement any unitary operation (and 
simulate any other circuit). A measurement on the output wires gives the result of the computation. Note that the 
QCM comes equipped with a canonical decomposition into a tensor product of small systems - the qubits (or qu-eíits 
in the general case of d dimensions). This is of course not only necessary to define a bàsic unit (like a bit for a classical 
computer), but also to decompose a circuit into bàsic gates. Both are indispensable for the uniformity of the circuit 
family: given a size of a problem (e.g. to factor a number given as a bitstring of length n) one should be able to 
efhciently (i.e. polynomial in n) find the necessary gate sequence. This would be impossible if the structure of the 
Hilbert space changed with inercasing n. The tensor structure ensures that this is not the case, and the QCM is 
thereby amenable to analysis in terms of the quantum complexity of an algorithm Q . 

We note here, as pointed out by many other authors jï] g], that all known physical models which can perform 
some form of quantum computation appear to be simulatable by the quantum circuit model. We therefore take this 
cmpirically verified hypothesis as a fundamental definition of what it means to build a quantum computer. Thus, for 
the issue of universality, some map between our physical system and the QCM must be established. In particular, 
one of the essential features of the QCM is a tensor product structure between different subsystems and therefore we 
shall impose this on our physical system. 



B. Universal Gate Sets 



In both the classical and quantum regime, universality results have been established that ensure the equivalence of 
ostcnsibly different models of computation which may employ different sets of gates. Studies of universality in the 
QCM have culminated in the definition of a universal gate set. A natural impulse is to make a close analogy to the 
classical case, i.e. to search for a finitc set of gates that exactly implements every quantum circuit efficiently. However, 
these attempts are doomed to failure: obviously a discrete set of gates can't be used to implement an arbitrary unitary 
operation exactly, since the set of unitary operations is continuous. Powerful results have nevertheless been established 
for the approximate simulation of quantum circuits with a set of finite gates, and we will now briefly review the steps 
leading there. A more detailed account can be found in JhJ. 

To quantify the aceuracy of an approximation, one needs a distance on matrices. On the finite N dimensional 
matrix space, any mètric is basically as good as any other (i.e. they differ by only constant factors). For example, 
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we can use the trace-norm d(U,V) = J 1 - -^Rc [Tr(UtV)], or the norm d(U, V) = max^\\(TJ - V)|*)|| (|*) is a 

normalized state). A matrix V is then said to approximate a transformation U to accuracy e if d(TJ, V) < e. 

The methodology employed in most universality proofs and constructions starts with a set of gates and then operates 
within the unitary group generated by their repeated applications to show that the whole of SU(N) is generated, 
where N is the dimension of the spacej 1 . Historically the first universal quantum gate sets were three-qubit gates 
(Deutsch ||). Earlier work by Toffoli jÏ2| on classical circuits has introduced a classical universal three-bit gate - the 
Toffoli-gate - which is reversible. Classical circuits based on the Toffoli-gate can simulate all other classical circuits 
efhciently and in a reversible way Deutsch generalized this gate to a unitary three-qubit quantum gate and proved 
its universality. In 1994-96 several groups improved on Deutsch's circuits. DiVincenzo []ïl"f demonstrated how to 
construct Deutsch's three qubit gate using only one and two qubit gates, and thus demonstrated that the physically 
unrealistic three qubit interactions of Deutsch are not necessary for universality. Barenco et al. |Ï3| proved that single 
qubit gates together with the controlled not (CNOT) allow one to exactly implement any unitary transformation. This 
is now so often taken as the paradigm for universal computation that we may refer to it as the "standard model" . 
Another memorable universality result is due independently to Deutsch, Barenco and Ekert JÏ4| and Lloyd fïj^ in 
1995. They prové using Lic-algebra arguments (see next section) that almost any two-qubit unitary gate is universal. 
Almost any is to be taken in the genèric sense: in the 16-dimensional manifold of possible two-qubit unitary gates 
acting between n qubits, the set of gates which do not yield a universal set of gates is a lower dimensional manifold. 
They implicitly assume that the genèric gates (given by Hamiltonians Hij) can be applied between all qubit pairs 
i, j and in both directions (i.e. as Hij and Hji), i.e. they assume that some sort of permutation (or exchange) 
gate is given. Howcver many of the accessible interactions in physical systems are not of this sort - nature is often 
non-generic\ 

To quantum computo in the prcscncc of (inevitable) errors, quantum information has to be protected against noise. 
To this end several approaches have been developed, among which are the theory of quantum error correcting codes 
(QECC) jÏ6| and the theory of decoherence-free subspaces and systems (DFS) [ p"7|fl8[ ]. In these approaches quantum 
information is encoded into parts of the system Hilbert space. The theory of fault-tolerant quantum computation 
shows how to achieve universality in the presence of errors and imperfect quantum gates. The bàsic idea is to work 
directly with encoded universal gates that act on an encoded space, whether this derives from a quantum error 
correcting code or from a decoherence-free subspace or subsystem. Considerations of fault-tolerance have led to a 
number of specific universal gate sets. These gate sets are useful because it is known how to implement them within 
the context of fault-tolerant quantum error correction. We will mention only one of them here ]Ï9|], namely the set 
consisting of {H, P, CNOT, n/8} where 
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For practical implementations, the degree of locality of the gates is important. For QECCs it has been shown how 
to implement such an encoded discrete universal set of gates using local gates on the physical qubits |f(|. These 
gates take encoded states to encoded states, but may take the state outside of the encoded subspace during the 
computation. In order to apply this approach, the speed of gate operations is critical: it has to be less than the 
error rate. In the theory of computation on DFSs, Bacon et al. jïj (for a specific class of DFSs) and Kempe et al. 
H (on a broader class of DFSs) have shown how compute on DFS with local Hamiltonians only. This work shows 
that the two-body exchange Hamiltonian alone is sufficient to implement universal computation on a decoherence-free 
subsystem. This demonstrates in particular, that there is an encoding of the quantum information into a subspace 
such that the exchange Hamiltonian - which is a priori not universal in the sense that it can perform any unitary 
evolution on all qubits - does now become universal. We have termed this "encoded universality" . Ideas about encoded 
universality in the context of physical realizations of a quantum computer have been presented in j^,||,^0|,|I| and also 
implicitly in |p2| in the context of dynamical noise suppression schemes (so called "bang-bang" control) . 



1 A notable exception is the proof by DiVincenzo |TÏ| 
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C. Lie Analysis of Intrinsic Physical Interactions 



Every quantum gate is generated by a Hamiltonian. It is very convenient to recast the universality rcsults into a 
Hamiltonian language. This allows us to use Lie-algebra analysis to approach universality constructions. Physically, 
the implcmentation of a unitary quantum gate will come about duc to the implcmentation of some controllablc 
Hamiltonian for a specific time. 

Instead of supposing that we are given a set of unitary gates, we thus rather assume that we have control over 
a set of Hamiltonians which we can apply to our system for a controllable amount of time. Of course, there will 
be inaccuracies in the application of a Hamiltonian for any given exact time. These inaccuracies, however, are 
exactly what fault-tolerant quantum computation is designed for. Consequently it is lcgitimate to proceed from the 
assumption that we can exactly implement a given Hamiltonian for a given amount of time. The aim is then to a) 
rcformulate the analogue of universal sets of unitary gates Hamiltonians formulation and b) to provide efficient finite 
time approximations that result in implementable sets of universal gate operations. 

We can now ask what new gates (and new Hamiltonians) we can obtain by operating these given ones. A general 
theoretical procedure to combine known gates to construct new gates is well known, deriving from the Baker-Hausdorff- 
Campbell operator expansion. Accordingly, we can combine two Hamiltonians to construct new Hamiltonians using 
the following properties: 



In practice we will truncate the infinite series to approximate the relevant operations. A detailed analysis of this is 
provided in |Ï5|,p|]. As a consequence of these composition laws scalar múltiples, sums and Lie-commutators i[A,B] 
can be obtained out of the given Hamiltonians) . Note that these operations correspond to closing the set of allowed 
Hamiltonians as a Lie-algebra. To approximate to a given accuracy, say exp(i(aA + /3B)), by a sequence of gates 
with generators A and B one just truncates Eq. (^) for large enough p. This allows us to approximate every unitary 
element in the Lie-group corresponding to the Lie-algebra generated by the given Hamiltonians using controllablc 
coupling constants. The inverse of this statement is also true, namely all gates that we can build of a set of given 
interactions are generated by elements in the Lie-algebra. 

It is now easy to define a universal set of generators H in this setting: H is universal if the Lie-algebra generated 
by the elements of H contains su(N), where N is the dimension of the underlying space and su(N) is the Lie-algebra 
of the special unitary group. This implies that the set of unitary gates obtainable via successive application of gates 
deriving from H is dense in the group SU (N) . 

This general approach provides a systematic way to construct universal sets of gates, but so far it says nothing 
about whether or not the resulting gates can be implemented efficiently. We remind the reader that we are using 
the term "efficiently" in the quantitative sense defined above, namely in terms of whether implcmentation can be 
made with only polynomial overhead in the number of discrete gates required. The Trotter expansion, Eq. (||), is 
guaranteed to be aceurate in the limit as the number of terms p becomes very large, but clearly a critical practical 
issue is how the truncation seales with the size of the quantum circuit, n. Thus in order to actually use this approach 
to universality for realistic physical implementation, it is crucial to have bounds on the length of the gate sequences 
approximating a certain gate in terms of the desired accuracy for a circuit of given size. This is all the more important 
if one universal set is to be replaced by any other with only polynomial overhead in the number of gates applied, for 
otherwise the complexity classes would not be robust under the exchange of one set for another. The Solovay-Kitaev 
theorem^ establishes the equivalence of universal sets, and provides bounds on the length of gate sequences for a 
desired accuracy of approximation. In short, this theorem states: 

Theorem (Solovay-Kitaev): — Given a set of gates that is dense in SU(2 k ) and closed under Hermitian 
conjugation, any gate U in SU(2 k ) can be approximated to an accuracy e with a sequence of poly [log(l/e)] 
gates from the set. 

It follows that any quantum circuit (performing any desired unitary operation) can be simulated to arbitrary accuracy 
efficiently by a circuit with gates from a universal set. In the context of a universal set of Hamiltonians H it implies 
that all the gates in the Lie-group can be obtained efficiently by an appropriate choice of the control parameter. This 
powerful result lays open the route to searching for radically different universal sets of operators. 



2 For details and proofs see W 
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III. FORMALISM 



Designs for quantum computers can come from a variety of undcrlying physical systems. A major challenge for 
futurc quantum engineers is to find the physical systems that allow for a universal set of gates. In all the physical 



implcmcntations proposed only particular Hamiltonians can be turned on and off. As mentioned in Sec. II B, genèric 
two-body interactions can suffice to implement any unitary transformation. Nature, however, is not genèric. Most 
interactions in real physical system possess symmetries that place these interactions squarely on a lower dimensional 
"non-generic" manifold where they are not universal. Thus, while allowing in principle for universality based entirely 
on two-body interactions, this approach does not provide a general prescription for universality given a specific set of 
intrinsic interactions. 

We propose a different approach here, namely to start from the natural interactions given by the physics of the pro- 
posed qubit system, and then to investigate the undcrlying potential of these interactions for "encoded universality" . 
The rationale behind this is that the given interactions - although per se not universal - become universal within a 
subspace. Encoding into this subspace allows to reduce the number of different Hamiltonians needed for universality, 
at the expense of spatial resources. This approach allows one to take the natural physical constraints of a specific 
implemcntation into account. For a given physical system, the natural interactions of choice can then be dctcrmincd 
by various factors, including their speed, the ease with which they can be implcmcntcd, any cnginccring or material 
constraints, and their robustness towards decoherence processes. 

Here we detail this new paradigm of "encoded universality" . We start by reviewing the relevant results from our 
previous work addressing quantum computation on decoherence-free subsystems and by summarizing the mathemat- 
ical framework. We illustrate the general theory with a summary of our previous results on the isotropic exchange 
interaction before presenting the new results for the anisotropic exchange interaction in Section IV. 



A. Representation Theoretic Framework 



The general formalism employed to find universal encodings is largely related to the Lie-algebra framework for 
noiseless or decoherence-free (DF-) subsystems We briefly review this framework here. 

In the context of DF- subsystems we are dealing with operators in the interaction Hamiltonian H/ = ^ Q S Q (g> 
B Q . The Hamiltonian H/ describes the system-bath interaction which ultimately causes decoherence and corruption 
of the information stored in the system. The S a (acting on the system) generate an interaction àlgebra S, which is 
the àlgebra generated by the set Si = {7, S%, S2, ■ ■ ■ } by linear combination and operator multiplication. Because the 
interaction Hamiltonian is Hcrmitian, this àlgebra is f-closed. The àlgebra S splits into irreducible components as a 
consequence of the following general theorem from the representation theory of f-closed algebras (see, e.g., [M): 

Theorem 1 Let S be a \-closed àlgebra of operators acting on a space Tís- Then the space is isomorphic to a direct 
sum 

Hs~J2 C nj <8> C dj (4) 

in such a way that in this representation S and the commutant S' of S are decomposable respectiu ely, as 

S = © I nj ® M(C d > ) S' SÉ © M(€ n ' ) ® I dj (5) 

JEJ JGJ 



Here M(C") means the set of all linear operators from C™ to itself, and J € J runs over all the irreducible representa- 
tions (of dimension dj and with degeneracy n,j) of S. The commutant S' of S is the space of all operators commuting 
with S. It is itself an associative àlgebra, i.e. an àlgebra closed under linear combination and multiplication. States 
encoded in C™ J are completely immune to the interaction with the bath, because the interaction àlgebra S acts only on 
its "co-factor" C dj . Thus C n J is a noiseless or decoherence-free subsystem where information is intrinsically stabilized 
against the effects of the interaction with the bath, with no need for corrective action. It is important to note that the 
Lie-algebra generated by the S a via linear combination and Lie-commutator is different from the interaction àlgebra 
S. This point will become important later. 

One of the most important physical models of decoherence, which allows for decoherence-free subspaces and subsys- 
tems, is the model of so called collective decoherence. In the case of collective decoherence on n qubits the interaction 
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àlgebra S is spanned by S Q = X)T=i °a' w ith a = x, y, z. This àlgebra is well known in the quantum mechanical liter- 
ature, it represents a (reducible) representation of the Lie-algebra of su(2). Its block-structure according to Theorem 
|ï| allows to find the decoherence-free encoding of the information. 

Yet this encoding alone leads only half the way towards decoherence-free quantum computation. In order to harness 
DF-subsystems as quantum computers it has to be shown that one can perform universal computation on DFSs in a 
fault-tolerant way using local gates on the physical qubits. Inspecting the statement of Theorem [ï] more closely it is 
not hard to see that operations that take DF-states to DF-states must be generated by Hamiltonians which lie in the 
commutant S' of the interaction àlgebra S. For universal computation within a DF-subsystem we need to identify 
local gates in this commutant (these will be the Hamiltonians we can control) and to show that the Lie-algebra 
generated by these Hamiltonians contains su{nj) on each cofactor C™ J . 

In our work on fault-tolerant universal computation on decoherence-free subsystems we showed that the 
exchange interaction alone generates the Lie-algebra su on every such DF-subsystem. This result was obtained for 
any set of n qubits, n > 3 and settled the question of universal fault-tolerant computation on decoherence free 
subsystems with local interactions. 

Let us now switch viewpoints from consideration of the subsystem to consideration of the interaction. Our result 
then implies that if we restrict the Hilbert space of the system to the space of the subsystem only, the exchange 
interaction alone is universal. In other words, the encoding into subsystems provided a way to make the exchange 
universal, meaning that the set of exchange interactions E = {Ey : 1 < i < j < n} on n qubits is a universal 
set of generators over a subspace according to the defmition in Sec. [II C| ! We point out that while this result was 
achieved specifically for the decoherence-free subsystems that are defined by resistance to collective decoherence, the 
decoherence-free properties of the subsystems are not required for the universality result. 

The scenario of encoded universality is in a sense the "inverse" problem to universal computation on DFSs. In 
contrast to the DFS situation, we are now presented with a set of given interactions H on the system Hilbert space 
TLs- The operations we can implement from this set generate a Lie-algebra C(H). The Lie-group generated by this Lie 
àlgebra (via "exponentiation" ) is a subgroup of the unitary group U and thus is a compact Lie group. Compact Lie 
groups, along with the Lie àlgebra which generates the Lie group, are completely reducible^] and hence the Lie-algebra 
splits into a direct sum of irreducible representations (irreps) 

£(H) Sí Cj{nj) ® I dj (6) 

JEJ 

over the state space 

Here Cj is the Jth irrep of £(H). This decomposition is reminiscent of that of the commutant S' in Eq. (|^), and 
so we use the same variables for the dimension nj of the Jth irrep of Cj and its degeneracy dj. (Note that these 
variables are interchanged for the interaction àlgebra S.) 

The problem we face in order to make the set H universal can now be formalized as follows: 

Given the Lie àlgebra C(H) = £j(nj) <g) Id, find a component C raj (g> <D dj of the state space over which Cj(nj) 

contains su(nj). 

If this can be satisfied, then we have achieved universality on a subsystem. 

The exchange-only universality theorem we proved in |2| shows that the Lie-algebra generated by the exchange 
interaction generates su(nj) on each component C" J . This implies that exchange lends itself for encoded universality. 
We can encode into any of these components (corresponding to subsystems in the DFS framework). The problem of 
the existence and structure of "encoded universality" with isotropic exchange is thus solved in the affirmative. 

Furthermore, we can immediatcly infer the coding efnciency of the "encoded universality" obtained from exchange. 
As shown, e.g., in gj the number of encoded qubits k over the number of physical qubits n approaches 1 as n — > oo. 
This means that for larger and larger n the amount of redundaney in physical qubits that are needed to encode 



approaches zero. However, as was outlined in Sec. II A, we do have to introduce a tensor product structure via what 



we will call conjoining (see Sec. III C), and this will place a constraint on the encoding efnciency. 



3 A representation of a Lie àlgebra is irreducible if the action of the representation on its vector space does not possess an 
invariant subspace. 
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The tensor product structure is important in order to map the QCM onto our subspace . The convenient way to 
proceed in construction of this structure is to use "clusters" of subspaces over which exchange is universal. These 
clusters may be dclimitcd either by some natural length scale in the physical system, or by an arbitrary cut-off. For 
instan ce the smallest possible block size consists of 3 physical qubits (to encode two states - a logical qubit) (see Sec. 
III D). The exchange interaction allows us to perform operations both within and between the blocks. The encoding 
cfhcicncy with such a tensor product structure is then limited by the block size to a value less than unity, no matter 
how many qubits are employed. Choosing a higher cut-off, i.e., a larger block size does increase the coding efficiency, 
but this will always remain at a value less than unity for finite block size. We will give morè details of the conjoining 
procedure in Sec. III C. 

We will briefly review some specific examples for isotropic exchange in Sec. III D. In Sec. IV we will show that a 
similar result also holds for the XY-interaction, i.e., for an anisotropic exchange interaction. 



B. Interaction Àlgebra 7^ Lie- Àlgebra 

The àlgebra representation Theorem |ï| has given rise to some misunderstanding, which we now clarify. Let us 
look at the case of collective decoherence on n qubits. The interaction àlgebra S is spanned by S Q = 5Zi=i °aj with 
a — x,y, z, via linear combination and multiplication. The Lie-algebra spanned by the S a - we will denote it by C(S) 
- represents a (reducible) representation of the Lie-algebra su(2). This is different from S, since it is obtained by 
linear combination and taking commutators, i.e. in general C(S) C S. However, for compact groups (and the S a form 
a compact group), it is true that when expressed in the same basis as S, the Lie-algebra C(S) splits into the same 
block-structure as that of Eq. (@). 

As we can easily infer from the DFS theorem ]25|,^| , if the encoded operations are to preserve the information in 
the encoded subsystems, the gate Hamiltonians should be in the commutant S' of S. It is exactly these operations in 
S' that act entirely within the subsystems, as can be seen from Eq. (Q). Yet it turns out that the commutant of C(S) 
and the commutant of S are the same; in both cases it consists of the set of elements that commute with all the S a . 

The good news is that instead of all S, we only need to study the structure of the irreps of su(2) as a Lie-algebra C(S). 
It is very well known in quantum-mechanics since it is precisely the Lie àlgebra deriving from angular momentum. 
This will give us the noiseless or decoherence-free subsystems, and also the commutant, S' . Note that the Lie àlgebra 
su(2) arises specifically from the DFSs resulting from the collective decoherence model. Other decoherence models 
will yield different DFSs that may be defined by different Lie algebras. 

As is well known from the representation theory of su(2) (see e.g. pq| ) the commutant of C(S) is related to the 
symmetric group (permutation group) S n |2q| . The natural representation of S n on n qubits is the set of operators 
that permutes the qubits, i.e. if n S S n , then n acts on basis states as 

7T : |íl)|i 2 ) • ■ ■ \in) ► Ntt(1))|V(2)) ••·|V(n)) ïfc <E {0, 1} (8) 

Clearly the elements of S n commute with the permutation- invariant S a . But they - together with the associative 
àlgebra they generate - also constitute the set of all elements that commute with them. In other words the commutant 
of C(S) (this is equal to S' , the commutant of S) is given by the àlgebra spanned by linear combinations of elements 
of S n in its natural representation, Eq. (^). The associative àlgebra spanned by the permutation group S n gives the 
set of all operators on a subsystem, as can be seen from Eq. (|B|). In particular, we can find a set of generators for 
universal computation on subsystems among its elements. This means that for universal fault-tolerant computation 
all we have to do is to implement the permutation-group. It is also known that the permutation-group S n is generated 
by transpositions r.y (which just permute two qubits). These transpositions can obviously be implemented by local 
interactions - just switch on the exchange Hamiltonian. One might be tempted to conclude that the question of 
universal computation with local gates is thus settled! 

Unfortunately this conclusion is false. The transpositions generate the permutation group S n via multiplication. 
But the composition laws for given operations do not generate produets of allowed Hamiltonians. The allowed ways to 
compose operations from a set of bàs ic H amiltonians give only linear combinations and commutators, i.e. they close 



the Lie-algebra of the bàsic set (Sec. II C ). In other words, if our bàsic interactions are the exchange (transposition) 
of two qubits, we can implement the Lie-algebra generated by them, but not the associative àlgebra. 

For universal computation with exchange we had therefore to show explicitly B that the Lie-algebra generated by 
the transpositions Tij generates the whole Lie-algebra of su(nj) on a subsystem <D nj . 

In order to clarify this it is useful to discuss an example. Suppose that we are given a system of n qubits upon 
which we can enact one of the three collective Hamiltonians 



Ca=5>Í- (9) 
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The Lie àlgebra A generated by these interactions is simply that of su(2): 

[C a , C/3] — ie Q/ 3 7 C 7 . (10) 
This Lie àlgebra will thus be reducible to a set of nj dimcnsional irreps of su(2): 

n/2 

A= Aj{nj)»I dj , (11) 

J=0(l/2) 

where Aj(nj) is the n,j = 2 J + 1 dimensional irrep of su(2) which appears with degeneracy dj in the decomposition. 
It is important to realize that if a set of Hamiltonians represents an nj dimensional irrep of su(2), then the operational 
power of these interactions is not more than the operational power of any other irrep of su(2). In particular, it is 
no more powerful than the two dimensional irrep of su(2). An nj dimensional represcntation of su(2) acts on a 
nj dimensional space, but does not enact su(nj) on this nj dimensional space. Yet if we look at the associative 
àlgebra spanned by the C a and its decomposition into a block structure according to Theorem [ï], we see that on each 
co-factor of each block this àlgebra generates every possible matrix, and in particular, all matrices in su(nj) (it is the 
same associative àlgebra that we studied in the case of collective decoherence) . The point is that the power of the 
associative àlgebra is nevertheless not accessible via the fundamental physical composition laws. 



C. Conjoining and the Tensor Product Nature of Computation 



The above summary of representation theory for Lie algebras has been made solely in terms of some abstract 
(i-dimensional Hilbert space. We now admit two tensor product structures: first, on our physical system, and second, 
on the encoded qubits from which encoded universality will be constructed. 

The first tensor product structure is simply that implied by our physical system and is usually forced on us by the 
locality of physics. For example, if we are using the spins of single electrons on a quantum dot, the tensor product 
is just the natural one of these spin-qubits. The nature of this tensor product structure is also manifest in the set of 
interactions which will be present in the real world, and is rcprcscntcd in the Hamiltonians in H at our disposal. This 
tensor product is not relevant to our discussion except in the context that we expect the physical interactions we are 
dealing with to be local with respect to it. 

The second tensor product is the important tensor product relevant to building a map between the given interactions 
and the QCM. As mentioned before within an subspace for universality there is a priori no tensor structure present. 
In order to build uniform quantum circuits for quantum computation we need to introduce a tensor structure at some 



point ( Sec. [I A). In particular we are sceking to use a certain set of interactions to simulate a QCM and therefore 



we require that there be a tensor product structure corresponding to the the tensor product structure in this model. 

In order to define a tensor product structure, one must be able to identify subsystems with a tensor product structure 
between them. The exact method for establishing where this tensor product structure comes into play is arbitrary. 
From a physical perspective one is motivated to search for small encodings, i.e., we seek to minimize the number of 
physical qubits used to encode a logical qubit (or qudit if d-dimensional). The tensor product structure will then be 
induced between these encoded qubits (qudits) by separating them into blocks (e.g., of nearest neighbors), with the 
dcfming property that single encoded-qudit operations are possible within each such block. However the cut-off is in 
principle largely arbitrary and so our ultimate choice will depend on the balance between a desired coding efficiency 
versus the complexity of gate sequences required for this. These blocks of encoded qubit (or qudit if d-dimensional) 
systems will form the bàsic factors of our ultimate tensor structure. We term this process "conjoining" the encoded 
qudits Within each block we can implement any encoded operation as a consequence of our universality proofs. 

A subtlety arises when we want to perform encoded operations between the blocks. This has to be addressed 
separately from the operations within a single block. In the case of the exchange interaction we have the property 
that the space formed by the two separate encoded qudits is itself part of an encoded higher dimensional system. 
This is a consequence of the structure of the subsystems, as shown explicitly in the proof in This means that the 
exchange interaction implcments any operation within this bigger subsystem, including in particular, any operation 
we want to perform between the blocks. A similar statement is true in the case of the XY-interaction. 

Finally, given a certain encoded universality construction, the issues of preparation and measurement on the encoded 
information must be established. Generally speaking, some method of efficiently generating states with a large overlap 
of the encoded information must be possible. Similarly, a method for extracting some bit of information from the 
encoded information must also be available. These are important issues which need to be addressed within the physical 
feasibility of a given encoded universality construction. 
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Finally we note that the formalism we have developed here works only for a quantum circuit where measurement is 
not used as a fundamental process for enacting unitary evolution. Recently it has been demonstrated that universal 
quantum computing can be achieved using only measurements on a particular prepared entangled state [p7| . Our 
formalism does not deal with measurements used to construct approximate unitary evolution. 



D. Example - Universality from Isotropic Exchange Hamiltonian 

We will illustrate the theory of encoded universality with an example of our previous results on the isotropic 
exchange interaction || . The most general form of the Hciscnberg exchange interaction between two spins is the fully 
anisotropic spin-spin coupling 

H y = jXaÍaí + J%oÍaÍ + J?oÍ4. (12) 



L· iJ ~ "ij u x u x ~ °ij u y u y ~ 

„11 ; n , _ _ 

»J V U ~~ l 3 ' 



In the isotropic limit the exchange couplings between a given ij pair are all equal, i.e., — = = Jij, so that 



we are then dealing with the Heisenberg Hamiltonian 

Hl·lei = ^ Jij~Eiij. (13) 

To show how this intrinsic coupling can lead to encoded universality, we analyze its action on three qubits, assuming 
for simplicity that the there is no dependence of the amplitude J on the qubit indices, i.e., = J. 

The Lie àlgebra Ce generated by Ey on three qubits allows us to implement the Hamiltonians in the set 
{E12, E23, E13, T = i([Ei2, E23]}. A better basis for the Lie àlgebra is given by the set of operators Ho = 
E12 + E 23 + E13, Hi = j^j (E13 - E 23 ), H 3 = j2 (-2E 12 + E 23 + E13), H 2 = i[H 3 , Hi]. We then find that 



for all a and 



[H o ,H Q ]=0 (14) 



[Ha, H/3] — ie a [3~í~H-j (15) 



with a, /3,7 € {1, 2, 3}. Ho is an abelian invariant subalgebra of this Lie àlgebra and thus factors out as a global phase. 
The set {Hi, H2, H3}, on the other hand, act as the generators of su(2). Thus the exchange interaction between three 
qubits can be used to implement a single encoded-qubit su(2). More precisely, we find that the H Q for a € {1,2,3} 
generates the àlgebra 

Cf = (£i®/ 4 )®(£2®/ 2 ) (16) 

where Cd is the d-dimensional irrep of su(2). Note the degeneracy of the corresponding irrcps. Corresponding to this 
decomposition the H Q act as 

H Q = 4 © Qa- a ® (17) 

for a € {1,2,3} where O4 is 4-dimensional zero operator (the 1 dimensional irreps all act as 0) and I2 is the two- 
dimensional identity operator. The action of the exchange is thus identical to that of an su(2) operator on a single 
qubit when working over the encoded space defined by the above decomposition. If we encode our logical qubits as 

|0 L > = -^(|010>-|100» 

|U) = |001) - y||010) - y||100) (18) 

then we find that the action of H Q is hcr a . Due to the degeneracy of the irrep, other encodings are also possible. Note 
that these states are nothing but J = 1/2, J z = +1/2 total angular momentum states of 3 spin-1/2 partides with a 
given projection along a certain axis, and can thus be found using elementary addition of angular momentum 

Having shown that the action of the exchange interaction on three qubits can produce the effect of a single 2- 
dimensional representation of su(2), we now conjoin the qubits, i.e. induce a tensor product structure between blocks 
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of three qubits in order to simulate a quantum circuit. We thus choose an encoding scheme in which a single qubit is 
identified with 3 physical qubits. This is not the only choice of tensor product structure. In fact, any tensor product 
between sets of k > 3 qubits can be used to construct a universal gate set B. Once one has introduced a tensor 
structure for the encoding, it is necessary to show that the natural couplings of the system can produce a non-trivial 
action between the tensor components of this encoded tensor product structure, i.e., show that we get universal gates 
on the encoded qubits. 

In the above case, where we have shown that we can obtain full control (sn(2)) over an encoded qubit, what we now 
need to show is that a non-trivial action between the encoded qubits can be enacted. This is nothing more than a map 
from the encoded universality to the fully universal set of physical gates consisting of single qubit gates supplemented 
by a non-trivial coupling between the qubits (like the CNOT, for example). This point of the universality proof is 
typically the most daunting, but really amounts to nothing more than understanding the Lie àlgebra generated by 
the interaction over two tensored encoded qubits. 

For the exchange interaction we have shown explicitly that nearest neighbor interactions can be used to produce 
such an action In particular, we find that the effect of the exchange interaction on 6 qubits becomes 

Cf = (£ 5 ® h) © (£9 © h) © (£ 5 © h) © (Ci © I 7 ) (19) 

When we conjoin two encoded systems, the Hilbert space of the full system contains a subspace spanned by the two 
encoded systems. It is the action on this conjoined subspace that we must then show can act nontrivially to produce 
an encoded two-qubit gate. When two of the above three-qubit codes are conjoined, we find that the resulting tensor 
product subspace lie entirely within the second and third irreps of the decomposition, Eq. ( |Ï9| ) . This in turn implies 
that non-trivial interactions between states in the encoded irreps are possible. In fact the original encoded single- 
qubit operations su(2) are also contained within this decomposition. The most important point, however, is that on 
our established tensor product between encoded qubits we can indeed couple these encoded qubits in such a way as 
to achieve a map to the unencoded, fully universal set of gates. Thus we have established that with the exchange 
interaction, this encoding can efficiently (to within a factor of three in spatial resources) simulate the unencoded set 
of universal gates. 

While the above example provides a systematic procedure to arrive at encoded universality, it is not unique. 
Indeed, there is a great deal of arbitrariness in how to produce a universal quantum computer using encoding. Our 
scheme of taking small encodings and then coupling them by conjoining is by no means the only method to obtain 
a universal quantum computer. An interesting alternative approach that exploits topological aspeets is also under 
active consideration, see e.g., [^8 29|. However, in dealing with real physical implementations, we expect that the use 



of small encodings will be very convenient. The principle of conjoining small encoded qubits therefore provides a very 
useful and practical route. 



E. Efïiciency of Encoded Gate Implementation 

We now address the question of how efficiently these encoded gates can be implemented with discrete operations. 
Consider the encoded su(2) defined by the àlgebra of Eq. (|Ï5|). The Solovay-Kitaev theorem telis us that any encoded 
single-qubit operation can be implemented efficiently with some sequence of discrete gates e íHt ™ , where t n is the set 
of discrete times over which the exchange operations are switched on. However, it does not provide explicit vàlues 
for t n . There does exist a brute force route to such sequences of discrete gates, namely via truncation of the Trotter 
expansion, Eq. (^). However there is absolutely no guarantee that this expansion will be efficient. 

In recent work Q we have shown that while there does not yet appear to be an analytic route to determination of 
optimally efficient discrete gate sequences, it is nevertheless possible to use numerical methods to search for efficient 
sequences. Whether the resulting solutions are optimal or not is not known, but in the example of isotropic exchange 
the numerical results are striking. In particular, the numerical search in Ref. [Q showed that there exists a sequence 
of 19 discrete exchange gates that can achieve a nontrivial two-qubit gate between encoded blocks (the encoded 
CNOT), while encoded one-qubit gates require no more than 4 elementary exchange interactions. Several numerical 



approaches to search for efficient discrete gate sequences have been suggested 130-34 . Methods relying on local 
invariants of quantum states allow determination of discrete gate sequences for two-qubit encoded gates up to a one- 
qubit rotation (3^] . This was the procedure employed in Ref. where two such invariants were combined to yield a 
function whose minimal solution corresponds to a discrete gate sequence. Clearly there are many other approaches to 
the determination of optimal discrete gate sequences that could be explored, including factors such as pulse shaping 
and genètic algorithms. 
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F. A Useful Criterion for Demonstrating Non-Universality 

The above discussion has focused on how one can find universal encodings deriving from a given physical interaction. 
While we cannot always guarantee that this is possible and in general will need to proceed on a case by case basis 
to find such encodings, formal conditions for the "inverse" problem can nevertheless be made. Thus, in certain cases 
we can indeed answer the question of whether a set of interactions cannot be used as a universal quantum computer, 
even in the presence of arbitrary encoding. We present here a useful criterion for testing this, allowing detection of 
Lie algebras which are not universal. 

Suppose one is given a set of Hamiltonians H„ which can be implemented in a Hamiltonian control sequence on n 
subsystems. Let C n denote the Lie àlgebra which can be generated by H n , and let g{n) denote the number of lincarly 
independent operators in C n . We call g{n) the subsystems growth function. 

Theorem 2 A growth function g(n) which is polynomial in n is not universal on a quantum circuit model. 

Sketch of the proof: The bàsic idea behind this theorem is the realization that a quantum circuit model on n subsystems 
has a state space which grows exponentially in n. Therefore performing unitary operators on this space is equivalent 
to generating elements of an exponentially growing Lie àlgebra. 

We will illustrate the above by giving a set of Hamiltonians which is not universal for any encoding of the information. 
Consider the set of gates generated by Hamiltonians in the set 

C = {<7« 4'V£ +1 >} (20) 

where o~a is the ath Pauli matrix operating on the ith qubit tensored with identity on all other qubits. We will now 
show that, even with the help of encoding, this set of Hamiltonians is not universal. 

Define the operation 1Ajk,a,p — o~a^ T\i=j+i where j < k and a, fi G {x,y}. We claim that the operators 

in the Lie àlgebra generated by O' are all linear combinations of the form M.jk, a ,p plus the single qubit at^ ■ Notice 
that this is truc for n = 2. We will prové the result by induction. First we note that because our generators are made 
up of Pauli operators, we need not worry about linear combinations of operators, but only need to worry about the 
operators which can be generated by commutation. Let C n denote the Lie àlgebra on n qubits generated by taking 
commutators in O'. For example £2 = {^-I2,x,x>^-i2,x,y!^-i2,y,x^^i2,y,y, } as claimed above. Assume that 

£„ = {M.jk <olt p, 1 < j < k < n, a, G {x,y}} U {cri l \ 1 < i < n}. First notice that taking commutators of elements 
of £„ and o~i^ only produces elements in £„: the only elements which cr^ do not commute with are M^q,.^ and 
Mfcj iQ , p and this commutation only serves to flip the value of a or 0. Finally, note that taking the commutator 

' (i) 

between elements of £„ and o~ x o x can only generate elements that are linear combinations of the ~M.jk,a./3 and 
o~i l \ To see this, first note that the only nontrivial commutators are those which occur with the operators, which 
just produce elements in £2. Further commutators between a x a x and M.jk, a ,p only create 'M.j'k' ,a' ,beta' which 
are one qubit larger or smaller. Thus we have proved that the Lie àlgebra generated by elements of O' is spanned by 
the set of linearly independent operators in £„ = {M.jk lCll p, 1 < j < k < n, a, G {x, y}} U {o~i l \ 1 < i < n}. 

Let us count the operators in C n . There are n operators and 4(^) ~M.jk, a ,f3 operators. Thus the growth function 
for this Lie àlgebra is g(ri) = n + 4 ^ "^" 2 1 · ) ^ = 2n 2 — n. This growth function is polynomial in n and thus via Theorem 
^| this set of operators is not universal. 



IV. ANISOTROPIC EXCHANGE 



Anisotropic Heisenberg spin couplings arise whenever there is some preferred dircction in space along which the 
coupling is stronger or weaker. This could be due to, e.g., asymmetries induced by donor àtoms in solid-state arrays of 
àtoms coupled via their nuclear spins |35 
of the spins while the coupling in x- and y-direction is equally strong: 

{HxYhs = ^K4 + °y v ) = JíúAíú- (21) 

This situation is relevant to the proposal of solid-state quantum computation using quantum dot spins and cavity 
QED §3. 



36 1. The XY-interaction arises when there is no coupling in the z-direction 
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We will now analyze the power of the XY-interaction as a universal gate via encoding. Several recent works 
have addressed ways in which the XY-interaction can be made universal under the addition of either single-qubit 
operations [ |37|j3l| , or single-qubit interactions (3^]. In contrast to these proposals, our focus here is on deriving 
universality without any additional single-qubit interactions or operations. Our presentation here is designed not only 
to illustratc the general formalism and methodology of proof, but also to make the approach to find encodings that 
provide universality with a single physical interaction, transparent and applicable to other types of interactions (see 
also pof ). Technical details of the proofs are relegated to the appendix [A. 



A. Lie Àlgebra for Anisotropic Exchange 

The set H = {A;j : 1 < i < j < n} generates the Lie-algebra £ To help us find the splitting into irreps Eq. (^|) let 
us identify (by inspection) some elements in its commutant £: 

n 

s, = 5>* x = ns =1 «4 (22) 

i=l 

We claim that the associative àlgebra Ai generated by these two elements via linear combination and multiplication 
is identical to the commutant £ of £ Note that the commutant is always closed under linear combination and 
multiplication and hence is an associative àlgebra. Together with the identity operator, Ai obviously constitutes a f 
closed àlgebra and hence satisfies the condition of Theorem [|. It splits into irreps according to Eq. (j^) . Assume for 
now that £ = Ai. Then the splitting of C and the splitting of Ai (as well as that of its commutant Ai') will be dual 
to each other: 

£ SÉ £j(nj) I dj M SÉ I nj M(C dj ) M' Sí M(C nj ) ® I dj (23) 
JeJ JeJ JeJ 

in the same basis. In particular the splitting of Ai would give us the encoding into subspaces over which H is 
possibly a universal set. Note that to show that H is indeed universal over these subspaces we still have to prové 
su{nj) Ç C,j(nj). 

Now consider that our claim of equality between £ and Ai was wrong, and that actually £ D Ai holds. This 
implies that C C Ai' . This would mean that the splitting of Ai' is coarser than the splitting of £, since adding 
elements to C just "joins" several previously separate factors (J) C nj (8 C into one bigger block C 7ljdj ® This in 
turn would mean that C in this basis would have a sub-block structure. In particular, £ could then not generate a 
full su(nj) in this block because it doesn't mix certain parts of the bigger block. 

Our strategy is then the following. We will first identify the splitting of Ai and its associated subspaces. Then, 
using an inductive proof, we shall show that over these subspaces C contains su(nj), thus allowing for universal 
computation. As a by product this will also prové our claim above that £ = Ai does indeed hold. 



B. The Àlgebra M 

Let us study the àlgebra generated by the two elements S z and X in Eq. ( p^ ) . They do not commute so we expect 
some irreducible representations of higher dimension. The element S z alone is already diagonal in the Standard basis, 
the eigenvalue of a bit-string with i ones and n — i zeros is n — li. The element X turns every zero into a one and 
vice versa. Let us rearrange all bit-strings into the following order 

B = {000 ... 0, 111 ... 1, 100 ... 0, 011 .. . 1, 010 ... 0, 101 ... 1, . . .} (24) 

i.e. listing elements with increasing number of ones and each followed by their transform under X. In this basis S z 
and X are block-diagonal with 2-by-2 blocks of the form [n— 2i)a z (for S z ) and a x (for X), where i = . . . [n/2\ is the 
number of ones in the corresponding first basis bit-string of the particular block. Obviously all blocks with the same 
i are the same, and there are (™) of these. By linear combination and multiplication we get the whole matrix àlgebra 
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M(C 2 ) on each blockg. In the representation of M (H) we have J = {0,1,2,... [n/2\}, dj = 2 and nj = (")^ 

The largest encoding - assuming we can prové universality of L on these spaces - will be of dimension nj — Q n _"y 2 ) 
for odd n and ( n / 2 -i) ^ or even n l3 ^ e dimensions of the first few encodings are given as follows: 



n 


1 


2 


3 


4 


5 


6 


7 


8 


max nj 


1 


1 


3 


4 


10 


15 


35 


56 



Thus, for rt = 3 we can encode one qutrit, for n — 4 we can encode two qubits, for n = 5 one can either encode a 
single ten-level system, or two five-level systems, or one can use only eight of the ten available states to encode just 
three qubits, etc. As n gets large, the rate of encoding approaches one (i.e. log 2 maxnj — > n), i.e. the encoding 
approaches unit efhciency. 



C. Example - One Qutrit Encoded into Three Qubits 



The smallest space into which we can encode morè than one state is the space of n = 3 qubits. This encodes a 
logical qutrit. The àlgebra M. and its commutant C C Ai' split as 

M = h <8> M(C 2 ) 8 h <8> M(C 2 ) M' = M(<tf) ® I 2 ® M(€ l ) ® I 2 - (25) 

We will show "by hand" that C generates all of su(3) over the second product. According to the previous section we 
can chose the basis for our encoded qutrit as 

|0 L > = |100) IU) = |010) \2 L ) = |001). (26) 

The explicit action of H on this encoded space is 

/o o o\ /00l\ 

(27) 

where is an operator on the encoded space. Taking commutators of two of these encoded operators gives us ia y 
on each two-dimensional subspace. Further commutators can easily be seen to generate all of su(3). Hence H is a 
universal set on the encoded qutrit. 






We can now introduce conjoining of the encoded qutrits (see Sec. [II C) in order to make use of them as bàsic 
units in an encoded circuit. As outlined above, for this to be useful, we have to show that we can implement some 
non-trivial coupling between the qutrits such that the full su(3 x 3) — su(9) can be generated between two conjoined 
blocks. For example, it can easily be calculated that for the 6 qubits made by conjoining the two 3-state codes, the 
operator [L4i6, A15], A12] produces a coupling which preserves the conjoined coding and non-trivially couples the two 
encoded qutrits. By putting together encoded qutrits (conjoining) in this manner, the XY model can therefore serve 
as a universal quantum computer. 

Note that in general gates between encoded blocks do not constitute a problem once we have proved universality 
of H for higher dimensions. The 9 basis states obtained by conjoining the 2 encoded qutrits give states with exactly 
2 ones (and 7 zeros). In the 6-qubit space they are all in the same subspace, namely the space generated by states 
with exactly 2 ones. They span a 9 dimensional subspace of this 15 dimensional space. The general proof (see below 
and Appendix) will show that it is possible to implement su(15) on this space, and in particular the sub-algebra 
su(9) ffi su(6) C sit(15), where su(9) acts on the 9 states obtained by conjoining two encoded qutrits. Hence we can 
obtain any encoded gate between them. 



4 There is one exception in the case of even n. Over all the ( n ™ 2 ) states with equal number of zeros and ones S z acts as 0. 
This means that over this subspace the àlgebra Ai is generated by X only, is thus Abelian and splits into one-dimensional 
irreps in the basis where X is diagonal. It is easy to see that this basis is formed by states |s) + X|s) (where X acts as 1) and 
\s) — X|s) (where X acts as —1), where \s) is a bit-string. Over this subspace M splits into two different 1-dimensional irreps 
with degeneracy ( n " 2 ) /2 as M = I( »^ /2 ® M(C X ) © 1^ ™j /2 ® M(C l ). 

5 and n n /2 = ( n / 2 ) f° r even n 
6 note: ( » )/2 < ( J_ ) forn > 1 
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D. General Proof 



We will now proceed to prové that for n qubits, where M! splits as[] 

M' = M(c("))®/ 2 (28) 

Je{0,l,...Ln/2j} 

the Lie-algebra C contains su(nj) = sm((")) on each component. Let us call each of the associated spaces S n (J). For 
instance S3 (1) is spanned by the three states in Eq. (p6[). In general we will make the canonical choice of taking our 
basis states to be those bit-strings which have #0 > #1^. Then for the anisotropic exchange, the label J actually 
counts the number of ones in the basis-bit-strings. 

We will proceed by induction, as in the corresponding proof for the exchange interaction to show that on each 
space H generates su(nj) independently on S n (nj), i.e. acting trivially on the rest of the space. We have seen in 
the previous section that in the case of n = 3 H generates an independent su(3) on S^l) (since all three generators 
annihilate 53(0) = |000) ) . Let us now understand how we build bigger spaces from smaller ones when we add a new 
qubit. To the basis states in 5„_i(J) we can either add or 1, which gives basis states in S n (J) and S n (J + 1) 
respectively. Turning this around this means that states in S n (J) come either from 5'„_i(J) - if a was added - or 
from 5„_i(J — 1) - if a 1 was added[]. 

For the inductive step let us assume that we can generate su(nj) independently on each subspace S n -i(J) of (n— 1) 
qubits. Within each space S n (J') of n qubits, we shall call those states that come from S n -i(J') by adding a zero, 
and 1 states that come from S n ~i(J' — 1) by adding a one. We call the states derived from going back two such steps 
00, 01, 10 and 11 states, by analogy to the TT, TB, BT and BB states used in the proof for isotropic exchange [||. 

Now let us analyze how a su on S n -i(J — 1) affects the spaces of n qubits^]. It propagates to a su on the 0-states 
of S n (J — 1) and the same su on the 1-states of S n (J). In the first step we are going to show how to eliminate the 
action of this su on one of these two spaces, keeping just the other one (independence step). This will give us an 
independent su on the states of each S n and another independent su on the 1 states. In a second step (mixing step) 
we will show how to "mix" these two su to obtain a full independent su(nj) on S n ( J). These two steps can be found 
in the appendix [a| which terminates the proof of universality of the XY-interaction on encoded spaces. 



E. Conjoining 



As we have seen on the example of the encoded qutrit in Sec. IV C it is possible to obtain all encoded operations 
between the blocks, which we have to introduce to obtain a tensor-structure. Assume we cut our blocks of size n 
qubits to encode a qudit. The XY-interaction generates any operation within a block. If we join two blocks, its basis 
states will give a subset of the basis-bit-strings in one of the encoded spaces over 2n qubits. (In fact, if we used the 
space S n (J) over the blocks then the bigger block will be a subspace of S2n(2J).) We can universally compute over 
this conjoined block and in particular obtain any gate in the tensor space of our encoded quditsf^] Hence we can in 
principle choose any cut-off we wish and increase our coding efficiency as much as we like. 

We now make some remarks on efficiency for the encoded XY. It is not clear how much overhead in time is required 
for actual implementation of this encoded universality with the XY-interaction. Numerical optimization studies 
similar to those we have made for the isotropic exchange interaction in pi and summarized in Section [II E have to 



be performed in order to find the minimal lengths of physical exchange gate sequences that can realizc specific gates 



7 the last component in the even case is M(c(™/ 2 ) /2 ) ® h ® M(C^" 2 ^ /2 ) ® h 

8 In the even n case we will have two spaces SÍ (n/2) and (n/2) of dimension ( n ™ 2 ) the first spanned by states of the 
form \s) + X|s) and the second by states of the form \s) — X|s). 

9 The only exception are the spaces S^(n/2) for even n. They can only be built by adding a 1 as the nth qubit to states in 
S n -\(n/2 — 1). This gives us a set of ( n ™ 2 ) strings |s). We then form the linear combinations |s) ± X|s) to obtain S„(n/2). 

We call this phenomenon "doubling" . Similarly, to build basis states in S„+i(n/2) from S„(n/2) we have to take all the ( n ™ 2 ) 
bit-strings with n/2 ones and add a 0. We call this procedure "rearranging" . "Doubling" and "rearranging" will require a more 
careful separate analysis in the inductive proof. 
10 By this we mean a particular operation in su(n.j) performed independently on this space. 

n If we used even n and one of the spaces S ± (n/2) for our encoded qudit, we have to apply a little more care with the 
argument. However the conjoined space will be formed by states in either of S2n(n) and hence there is no problem. 
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between encoded blocks. We note that several authors have recently shown that in a linear array, the XY-interaction 
is not universal when restricted to act between ne arest neighbors |40 - 43 only. In fact this result can easily be derived 
using the criterion for non-universality of Sec. III F by explicitly computing all possible commutators of nearest 
neighbor XY-exchanges. This set turns out to be polynomial in the numbcr of qubits. This result implies that the 
geometry of allowed interactions for our encoded universality result with the XY-interaction will be important in 
construction of numerical gate sequences. 



V. ERROR CORRECTION AND LEAKAGE 

We give here a brief discussion of the issues involved in making encoded universal computation fault tolerant. We 
do not pretend to give a complete treatment of this subject, but we provide some observations here that indicate that 
there is no fundamental obstacle to making our coded computation fully fault tolerant. 

Once we have a physical cmbodiment of a logical qubit (e.g., the sets of three spins with the states given in Eq. 
(|Ï8|)), and a procedure for implcmcnting CNOT and one-qubit gates on these qubits (J], then all the fault tolerant 
procedures developed recently for the Standard gate model for quantum computation can be applied |2íJ. The DFS 
coding would be concatenated within a stabilizer code, the appropriate transversal implementation of coded gates 
would be used, and appropriate error detection and correction would be implemented that limits uncontrolled error 
propagation. All of this carries over unchanged in the DFS setting. 

However, there is one aspect of error correction for which the DFS setting is quite different, and that is in the 
issue of leakage. Leakage is possible because the Hilbert spaces employed in quantum computation are larger than 
is assumcd in the ideal treatment. Ideally, the state of a qubit lives in a two-dimcnsional Hilbert space spanned by 
some states |0_l) and and the complete Hilbert space is just a tensor product of such Hilbert spaces. But for 

many physical embodiments of a qubit, the actual dimension of its Hilbert space is larger, even though we may have 
effective ways of restricting it approximately to just two of these dimensions. For example, an atom may have two 
low-lying energy levels that we identify as the Oi and lz, states; but of course the atom will have other excited states 
that we might label |2^) , |3i), etc. Unintended interactions and imperfections will surely, with some probability, 
bring the state of the qubit into these extra states, although we hope that the probability of this oceurring is small. 
This straying into extra dimensions is what is referred to as leakage. 

Leakage is actually a fairly likely error in the DFS schemes. For example, for the three-spin encoding of Eq. (|ï§|), 
The three spins naturally have an 8-dimensional Hilbert space, of which we are singling out two dimensions with 
quantum numbers J = 1/2, J z = +1/2. But very simple kinds of unintended interactions, such as a magnètic field 
affecting just one of the three spins, will cause the qubit to "leak" into the other six dimensions (those with quantum 
numbers J = 1/2, J z = —1/2, or with J = 3/2). So it is important to prescribe a scheme for detecting and correcting 
leakage to have a fully fault tolerant scheme. 



Data 



Ancilla I0> 




Measure 



FIG. 1. Circuit suggested by Preskill for leakage detection. 



Preskill [£9| has indicated a simple procedure for dealing with leakage. He suggests that, from time to time, each 
qubit be run through the simple circuit of Fig. The truth table of this circuit is to be 



\0l)\0l) ^ \0l)\Il) 

|li)|0i) "> |0£)|li) 

|2i)|0£) -> |2£)|0i) 
|leak)|0 L ) -> |leak)|0i 



(29) 



At the end of the circuit, the ancilla qubit is to be measured. If it is found to be a "1", the data qubit is in a vàlid 
state and it is unchanged by the circuit. But if it is found to be a "0", the data qubit is known to be in a "leaked" 
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state. In this case, the data qubit is replaced by a fresh |0l) state. This is usually incorrect, but it is an ordinary bit 
error, and is dealt with by the error-correction machinery at a higher level. 

The general strategy introduced by Preskill can still be used for the DFS-coded qubit, but the detailed procedure 
that he introduced will not work. From this point on wc will specialize our remarks to the case of coded computation 
using isotropic Heisenberg exchange. The problem is that the circuit given in Fig. ^ only has the truth table given in 
Eq. ( ^9|) given a certain assumption about how the quantum gates in the circuit work. In particular, it is assumed, 
for both the NOT and the CNOT, that if one of the inputs is leaky, the gate leaves the state unchanged: 



assumed CNOT action in Fig. 1: 
|0l)|0 l ) -> |0 l )|0l> 
\1l)\0 l )^\1 l )\1l) 
|2l)|0 l ) - |2 L )|0 L > 
|3 L )|0i) - |3i)|0 L ) 



(30) 



We find that it is not possible to implement a coded CNOT using the Heisenberg exchange that has this extended 
truth table, as it would require the gate to change the total spin quantum number of the state. Therefore, we must 
search for a different approach to leakage detection and correction. 



Data 



Ancilla I0> 



SWAP 
If 

Leaked 



Discard 



FIG. 2. Another leakage repair procedure better adapted to coded computation. 



Fortunately, we find that a modified leakage-correction circuit specification is permitted by the constraints of the 
exchange interaction; this circuit can then be used in the same way as the Preskill leakage-correction circuit within 
a complete fault-tolerant procedure. In fact, the circuit action we achieve, indicated symbolically in Fig. ||, can 
be used in a simpler way in fault tolerant quantum computation: the leakage is corrected without the need for any 
measurement. What we mean by "SWAP If Leaked" (SIL) in Fig. |is the following truth table: 

"SWAP If Leaked", Fig. 2: 
|0l>|0£> - |0£>|0 £ ) 

|U>|0 £ > -> |1 £ >|0 L > ^ 
\2 L )\0 L ) - |0 L )|2 L ) ^ 
|3l>|0 £ > -> |0£>|3 £ ) 

We note that if the data qubit is not leaked, it is not disturbed; if it is leaked, it is always replaced by \0l), permitting 
subsequent correction. (Of course, there is now no definite signal of when correction must be done. But if the rate of 
leakage is known, then error correction can be done with the appropriate frequency.) 

Actually, the function of "SWAP If Leaked" is achieved with a more relaxed specification: 

modified "swap if leaked" : 
|0l>|0z> - |0z)|0 L ) 

|1l>|0z>-|1z>|0 £ ) [ÓZ) 

|leak)|0i) — * |anything in 0l-1l space)|A). 

Here \X) is any state vector consistent with unitarity. This generalization will be needed in our DFS analysis. 

Now, we consider the implementation of SIL using the three-spin DFS coding discussed earlier. For this we adopt 
a specific notation for these states that will make the discussion easier: 

\N,n,J,J z ). (33) 
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Here N is the number of spins in the DFS (which is not really a quantum number, just a useful label), J is the total 
angular momentum, J z is its projection on the z axis, and n is a "principal" quantum number that labels the different 
states in the DFS (e.g., n runs from to 8 for the nine triplet states with N = 6, J = 1). 
In this notation our two coded qubit states are 



\0 L ) = |3,0,l/2, +1/2), \l L ) = |3,l,l/2, +1/2) 



(34) 



We will assume that a three-spin ancilla is available that has been freshly prepared in the coded state. Thus, we 
need to consider the operation of the SIL circuit on the eight basis states (the ancilla state is second): 



|3,0,l/2, +1/2) 
|3,l,l/2, +1/2) 
13,0,1/2,-1/2) 
13,1,1/2,-1/2) 
|3,0,3/2, +3/2) 
|3,0,3/2, +1/2) 
13,0,3/2,-1/2) 
13,0,3/2,-3/2) 



13,0,1/2, 
13,0,1/2, 
13,0,1/2, 
13,0,1/2, 
13,0,1/2, 
13,0,1/2, 
13,0,1/2, 
13,0,1/2, 



-1/2) 
-1/2) 
-1/2) 
-1/2) 
-1/2) 
-1/2) 
-1/2) 
-1/2). 



(35) 



3. through 8. are the "leaked" cases. To impose the constraint that the SIL circuit be implemented only with 
isotropic Heisenberg interactions, we must write these eight states in the basis of total angular momentum states for 
the six-spin block: 

1. |3, 0, 1/2, +1/2)13, 0, 1/2, +1/2) = |6, 0, 1, 1) 

2. |3,l,l/2, +1/2)|3, 0,1/2, +1/2)- |6, 1,1,1) 

3. |3, 0, 1/2, —1/2)13, 0, 1/2, +1/2) = l/\/2|6, 0, 0, 0) + l/\/2|6, 0, 1, 0) 

4. |3, 1, 1/2, — 1/2>]3, 0, 1/2, +1/2) = l/\/2|6, 1, 0, 0) + l/\/2|6, 1, 1, 0) 

5. |3, 0,3/2, +3/2)|3,0, 1/2, +1/2) = |6,0,2,2) 

6. |3, 0, 3/2, +1/2)|3, 0, 1/2, +1/2) = o|6, 4, 1, 1) + ò|6, 0, 2, 1) 

7. |3, 0, 3/2,-l/2)|3, 0, 3/2, +1/2) = c|6, 4, 1, 0) + d|6, 0, 2, 0) 

8. |3, 0, 3/2, -3/2) |3, 0, 3/2, +1/2) = e|6, 4, 1, -1) + /|6, 0, 2, -1). 

a-/are Clebsch-Gordan coeficients whose vàlues we will not need. Recali that in this basis, the Heisenberg interactions 
cannot change the J quantum number, and the action of the circuit must be independent of the J z quantum number. 
The "no-leak" cases 1. and 2. require, from Eq. (j3Ï|): 



(36) 



SIL\6, 0,1,1) = |6,0,1,1) 
SIL\6, 1,1,1) = |6,1,1,1). 

But because of the J z invariance of any circuit SIL made up only of exchanges, this implics 



SJL|6,0,1,S Z ) 
SIL\6,1,1,S Z ) 



|6,0,1,S Z ) 
|6,1,1,5,). 



This has implications for cases 3. and 4. We can use the decompositions 



|6,0,1,0) = 1/V2|3, 0,1/2, -1/2)13,0, 1/2, +1/2)- 
|6, 1, 1, 0) = 1/V2|3, 1, 1/2, -1/2)|3, 0, 1/2, +1/2) ■ 



1/V2|3, 0,1/2, - 
1/V2|3, 1,1/2, - 



-l/2)|3,0,l/2,-l/2) 
-l/2)|3,0,l/2,-l/2) 



(37) 



(38) 



(39) 



to write these out as 



SIL\3, 0,1/2, -l/2)|3,0,l/2, +1/2) = 
(1/V2)SIL(\6, 0, 0, 0)) + 1/2|3, 0, 1/2, -1/2) |3, 0, 1/2, 
SIL\3, 1, 1/2, -1/2)|3, 0, 1/2, +1/2) = 
(1/V2)5/L(|6, 1, 0, 0)) + 1/2|3, 1, 1/2, -1/2)|3, 0, 1/2, 



-1/2) + 1/2|3, 0, 1/2, +1/2)|3, 0, 1/2, -1/2) 
-1/2) + 1/2(3, 1, 1/2, +1/2)|3, 0, 1/2, -1/2). 



(40) 



We see that if we make the choicc 



SIL\6, 0,0,0) 
57X16,1,0,0) 



|6, 0,0,0) 
16,1,0,0) 



(41) 



17 



(note the sign change), then using these and the relations 



|6, 0, 0, 0) = l/y/2% 0, 1/2, -1/2)|3, 0, 1/2, +1/2) - 1/V2|3, 0, 1/2, +1/2) |3, 0, 1/2, -1/2) 

|6, 1, 0, 0) = 1/ v/2|3, 1, 1/2, -1/2) |3, 0, 1/2, +1/2) - 1/V2|3, 1, 1/2, +1/2) |3, 0, 1/2, -1/2), 1 ' 

we can plug back in to get 

SIL\3, 0, 1/2, -1/2)|3, 0, 1/2, +1/2) = |3, 0, 1/2, +1/2)|3, 0, 1/2, -1/2) 

SIL\3, 1, 1/2, -1/2)|3, 0, 1/2, +1/2) = |3, 1, 1/2, +1/2)|3, 0, 1/2, -1/2). ^ 



Here we achieve the generalized "SWAP If Leaked" requircment of Eq. (|3^ 

The last four cases are easier, because they are not constrained by the earlier choices. To deal with these, we 
introduce some new state expansions: 

5. |3, 0, 1/2, +1/2)|3, 0, 3/2, +3/2) = |6, 2, 2, 2) 

6. |3, 0, 1/2, +1/2)|3, 0,3/2, +1/2) = o|6, 7, l,l) + ò|6, 2, 2,1) 

7. |3,0,l/2,+l/2)|3,0,3/2,-l/2) = c|6,7,l,0)+d|6,2,2,0) 1 ' 

8. |3, 0, 1/2, +1/2)|3, 0, 3/2, -3/2) - e|6, 7, 1,-1) + /|6, 2, 2, -1). 

Note that the Clebsch-Gordan coefhcients a-f are the same as above. Then, we can fix the action of SIL to be 



S7£|6,0,2,7 Z ) = 16,2,2,7,) 
S7£|6,4,1,7 Z ) = 16,7,1,7,), 



(45) 



which will give us 



SIL\3, 0, 3/2, +3/2) |3, 0, 1/2, +1/2) = |3, 0, 1/2, +1/2) |3, 0, 3/2, +3/2) 
SIL\3, 0, 3/2, +1/2) |3, 0, 1/2, +1/2) = |3, 0, 1/2, +1/2) |3, 0, 3/2, +1/2) 

SIL\3, 0, 3/2, -1/2) |3, 0, 1/2, +1/2) = |3, 0, 1/2, +1/2) |3, 0, 3/2, -1/2) ^ 
SIL\3, 0, 3/2, -3/2)|3, 0, 1/2, +1/2) = |3, 0, 1/2, +1/2) |3, 0, 3/2, -3/2), 



as desired. This completes the analysis. We see that the specifications of Eqs. ( p8|]4l|j45| ) for the desired circuit are 
for a unitary transformation that is obtainable by only Heisenberg interactions, by the coded universality theorem of 
Q . We have not determined an explícit circuit for SIL; it might be interesting to use the techniques of Q to find one. 
Note that the specifications Eqs. ( p8| , [4Ï 45 ) for SIL only constrain its action on 16 of the 64 dimensions of the Hilbert 



space; its action on the remaining 48 dimensions is completely arbitrary, except for the constraints of unitarity. 

VI. OUTLOOK AND OPEN QÜESTIONS 

It is unclear at this point which general kinds of physical interactions are amenable to "encoded universality" . 
As we have seen, both the exchange and the XY-interaction give rise to a Lie-algebra that splits into degenerate 
irreducible representations. The degeneracy of this decomposition is the key to finding a suitable encoding. In 
physics, mathematical degeneracies are often linked to symmetries of a system, and indecd both the exchange and the 
XY-interaction possess a significant amount of symmctry. (The isotropic exchange is invariant under permutation of 
the Pauli matrices, and the XY-interaction is invariant under permuting a x and a y .) On the other hand, we know 



from previous work of Lloyd and others |15| that each genèric (and in particular non-symmetric) two-qubit interaction 
is universal without encoding if we allow bits to be flipped (i.e., the exchange gate). Note that both the exchange gate 
and the XY-interaction are invariant under the permutation of the specific two qubits on which they act. It would be 
interesting to know what the "encoded power" of other interactions not supplemented by an exchange interaction is, 
and also how this power is related to their symmetry properties. 

Another interesting question is related to the topology and arrangement of the physical qubits of the system. In 
many cases only nearest-neighbor interactions can be enacted. In our formalism we have allowed for application 
of the gate between arbitrary pairs of qubits. Thus, it would be important to know whether a restricted set of 
gates, for instance, only between qubits i and i + 1, has the same universality power. Clearly, in the case of the 
exchange interaction, we can limit ourselves to only nearest neighbor interactions, because we can decompose all 
permutations into such pairs. Thus exchange can be used to flip qubits at our convenience. However, as indicated 



above (Section IV E| ), this is not truc for the case of the XY-interaction. Whether for example certain two-dimensional 



layouts might be sufficient is still unclear. 
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A third open question is concerned with the overhead in time, i.e., the number of gates involving physical interactions 
switche d on for discrete time intervals that are needed to realize an "encoded" Standard gate. As discussed in 
Section [II E , there is currently no analytic route to finding the optimal gate sequences and we cannot even estimate 
the length of this in the general case. Establishing optimal discrete sequences is currently dependent on numerical 
searches in many dimensions. It would be highly desirable to develop systematic bounds on the length of sequences, 
even if an analytical route to optimization remains elusive. 
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APPENDIX A: DETAILS OF THE UNIVERS ALITY PROOF FOR THE XY INTERACTION 

Claim: Independence - Using operations in the independent su on S n —± (J— 1) and A„_i in we can generate 
an independent su on the 1-states of S n (J). 

Proof: Use two states in S n -i{J — 1) and a a z on them. This propagates to a a z on the 01-states of S n (J) and to 
a o~ z on the 00-states of S n (J — 1). Commute this o~ z with A n _i )n . Since A„_i.„ annihilates 00-states the net action 
on S n (J — 1) will be zero, whereas the net action on S n (J) will be ia y between one of the 01-states and another 10 
state. Commute again with A„_i „ to obtain an independent a z on the original two 01-states in S n (J). This new 
o~ z , now entircly acting only on S n (J), can be used via commutation with operations in the su to generate the now 
independent full su on the 1-states of S n (J) with the same arguments as in [|2j (using the Mixing Lemma of App. C 
there). 

When does this method not work? Only if there are less than 2 0-states in 5„_i(J — 1). If n— 1 > 3 this is never 
the case except for 5 Tl _ 1 (0), which just contains the |00 . . . 0) state. In this case S n (l) just contains one 1-state and 
sm(1) = is trivially implemented on itp|. 

Of course given independent su on the 1-states of S n we can immediately also get an independent su on the 0-states 
by subtracting the former from the coupled action of a su on a S n —i. 

Claim: Mixing - Given an independent su on each of the 0-states resp. 1-states of S n closing the Lie-algebra 
with A n —i tTl gives the full su on S n (independent). 

Proof: To get a ia y between a 0-state and a 1 state take a a z between two 10 states (or between a 00 and a 
10-state) and commute it with A n _x,n- To get a a z between these two states, commute with A n _i lTi again. This gives 
a su(2) connecting 0- and 1-states. The result follows from the Enlarging Lemma in [0 App. C. This proof works for 
all sizes of S n . 



12 Let us treat the case of even n here: "Doubling" - If we can implement su on S n -i(n/2 — 1) by construction this gives us 
coupled simultaneous su's on the two spaces S^{n/2) together with the same su on the 0-states of S n (n/2 — 1). The action 
on these latter states can be eliminated as before and we are left with these two coupled full su's on 5* (n/2). We will not 
attempt to eliminate the action of one of them, for practical purposes, if we need to use these spaces, we will only be able to 
use one of them. "Rearranging" - The problemàtic space is S n +i(n/2), more precisely its 0-states. A priory we only obtain 
two coupled smaller su's on this 0-space, and not the full su. However since rearranging provides a change of basis we do not 
have a block-structure and it is easy to see that commuting with A n -i, n mixes the two su's to give the full su on the 0-states 
of S„ +1 (n/2). 
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